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Chiral low-energy constants from lattice QCD 



1. Introduction 

The dynamics of QCD at low momenta can be described in terms of an effective theory, which 
encodes the spontaneous breaking of chiral symmetry 

SU{N f ) L x SU(N f ) R -» SU(N f ) v . 

Nf is the number of light quark flavours; we will consider both cases Nf = 2, 3. The effective 
Lagrangian in Euclidean space can be expanded in powers of momenta as [[l]-||] 

J2f z =JSfi 2) +JSf « + (1.1) 

where 

£>P = — Tr [dptfdpU] - -Tr \JHV + U^ f ] , (1.2) 

U G SU(A^) represents the pseudo Nambu-Goldstone bosons degrees of freedom, and ^# is the 
(Nf x iV/) quark mass matrix. The low-energy dynamics is parametrised by the so-called Low- 
Energy-Couplings (LECs) of the effective theory, which are not determined by symmetries. At 
leading order, the chiral Lagrangian involves two LECs, namely the pseudoscalar decay constant F 
and the chiral condensate £ (or equivalently B = L/F 2 ) l . At next-to-leading order, Jz^ 4 ' contains 
10 terms for the case Nf = 2 and 12 for Nf = 3, with corresponding couplings 

{Q -> / ; =i.. 7 A=i..3 (N f = 2), (1.4) 
{Q ^ L i=l . A0 ,H i=L . 2 (Nf = 3). (1.5) 

Once the LECs are fixed, chiral perturbation theory becomes a predictive framework and hence a 
powerful tool for understanding the phenomenology of QCD at low dynamics. A review of recent 
applications has been presented at this conference by G. Ecker 0. Those couplings must be ideally 
computed from"first principles": from this point of view lattice QCD is a promising approach, since 
-once the physical point has been reached (i.e. continuum, infinite volume limit and physical light 
quark masses)-it does not introduce any model dependence. 

For some decades lattice simulations have been performed far from the chiral limit and intro- 
ducing simplifications such as quenching. Chiral effective theory has been then extensively used 
to understand the results of lattice QCD, in particular to guide extrapolations to the chiral limit and 
to estimate the volume dependence of physical observables computed on the lattice. In the past 
years lattice computations experienced an important progress, due to an interplay of new theoret- 
ical developments, algorithmic improvements and increasing powerful computing resources. The 
effect of this progress is that now lattice unquenched simulations with Nf = 2 and Nf = 2 + 1 light 
dynamical flavours are approaching domains where a reliable matching with the chiral effective 
theory can be performed. Reliable in this case means that light quark masses are as close as possi- 
ble to the physical value, and that different sources of systematic errors can be at least in principle 
taken under control, namely finite-volume effects, lattice artifacts and renormalisation. 

1 In the following F, E and B will refer to the case Nf = 2; for Nf = 3 the notation Fq, Eq, Bq will be used. 



2 



Chiral low-energy constants from lattice QCD 



2. Low-Energy couplings from quark-mass dependence of pseudoscalar masses and 
decay constants 

The LECs can be extracted for instance by studying the quark-mass dependence of given 
observables in lattice QCD and match it with the predictions of the chiral effective theory. In 
this section I will summarize some of these predictions, in particular for the pseudoscalar decay 
constants and masses. The up and down quarks are clearly sufficiently light to formulate a SU(2) 
x SU(2) chiral effective theory where pions are treated as pseudo Nambu-Goldstone bosons and 
Kaons and etas are integrated out. In this theory, the LECS will depend on the higher energy 
scales m s ,m c ,mb, Aqcd- Alternatively, one can decide to treat also the strange quark as light and 
formulate a SU(3) x SU(3) effective theory, where also Kaons and etas are treated as pseudo 
Nambu-Goldstone bosons and the LECs will not depend on m s anymore. 

One of the simplest but deepest predictions which can be formulated at low energy is the 
so-called Gell-Mann-Oakes-Renner relation ph, which states that the square of the pion mass is 
linearly proportional to the light quark mass: 

M 1 n =M 1 = 2mB, (2.1) 

where m = (m u +m < i)/2. This corresponds to the leading order prediction in the effective theory. 
Higher order corrections can be computed systematically. In particular for the SU(2) x SU(2) 
theory one finds the NLO expressions [Q] 

m «= m+ 3^Hw> (2 - 2) 

M 2 ( M 2 



where A3 .4 are related to the scale-independent NLO LECs l^^ by the relation 

'3,4=ln(§A . (2.4) 

V / M= 139.6 MeV 

By matching the quark-mass dependence of F n and M % with lattice results, it is hence possible to 
extract the LO LECs F and £ (or equivalently B) and the NLO couplings I3J4. 
In this theory there is the possibility to treat Kaons as external matter fields coupled to pions, 
in analogy to what is done in the heavy-light meson chiral perturbation theory. The quark mass 
dependence of Kaon quantities can then be worked out explicitly in this framework [0]. 

In the SU(3) x SU(3) theory, the quark mass dependence of pion and Kaon masses and decay 
constants is given by ^ 

M 2 - = 2mB Q 1 1 + \i n - + 2mK 3 +K4V (2.5) 

M 2 = (m + m s )B j 1 + ^ + (m + m s )K 3 + K 4 X , (2.6) 
F n = F {1- 2\l n -Hk + 2mK 6 +K 7 } , (2.7) 
F K = F Q j 1 - Ifix - ^ K - 1^ + (m + m s )K 6 + K 7 X , (2.8) 
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Ml /M| 



with 

The Ki contain the NLO LECs: 

K 3 = —j (2L 8 - L 5 ) , = (m„ • m d • m s ) — y- (2L 6 - L 4 ) , (2.9) 
fir fb 

^6 = ^(L 5 ), K-j = (m u +md + m. v )^(L 4 ). (2.10) 

In this case, the matching with lattice results allows the extraction of Fq,'Lq{Bq),L 4 ,Ls : (2Lg — 
L5), (2L6 — L4). Moreover, the m^-dependence of the SU(2) x SU(2) LECs can be explicitly 
worked out [^] 

F = F ^l + S^L 4 -JI K + O(m^y (2.11) 

I = I 1 1 + 32^f U ~ 2Ji K — /ijj + OK 2 ) I , (2. 12) 

where M# is the Kaon mass in the limit m u ,md — > 0, m\ = m s Bo. In the large-A^ c limit one 
finds L/Lq,F /Fq,B/Bo — > 1, which corresponds to the so-called Okubo-Zweig-Iizuka (OZI) rule. 
Paramagnetic inequalities can be invoked to predict the sign of the deviations from the OZI rule, 
(F/F O -1)>O,(E/E O -1)>O[0]. 

Lattice simulations are performed on finite boxes with volume V = L 3 T, and finite- volume 
effects may play a relevant role. An important point is that chiral perturbation theory can be used 
also to infer information about finite-size scaling of physical observables [^|-p0||. In the asymptotic 
region M n L 3> 1, the dominant effect is of the form e~ M " L . The finite-size effects for the pion mass 
have been computed at two loops flnj] ; alternatively, resummations of asymptotic Luscher formulas 



[12] for pseudoscalar masses and decay constants have been investigated [13]. 

Lattice calculations are often performed in so-called partially quenched setup, where m sea ^ 
nival- It is possible to extend the effective theory formalism to this case and obtain predictions from 
partially quenched chiral perturbation theory [ p"4| ] which can be matched with the lattice results. 

Finally, lattice simulations are performed at small but finite lattice spacing a; ideally one 
should first perform a continuum extrapolation of lattice results and then match with the chiral 
effective theory. In alternative, one can formulate the chiral effective theory at finite lattice spacing 



by adopting the Symanzik formalism [ ]15[ , |lq ]. The price to pay is the appearing of extra couplings 
in the effective theory which have to be determined by the matching. Notice that for particular 
discretizations, where flavour symmetry is broken at finite lattice spacing, for example for the so- 
called staggered fermions with fourth root prescription, the limits a — > and m — > have to be 
taken simultaneously. Lattice results must be then matched with the so-called staggered chiral 
perturbation theory Jl7|]. 

In several recent lattice studies with Nf = 2 dynamical fermions the quark mass dependence of 
M K and F K has been matched with SU(2) x SU(2) chiral perturbation theory in order to extract F, £, 
1$, I4. The relevant parameters of these simulations are summarised in Table |l|: the discretization 
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Collaboration 


Dirac operator 


a (fm) 


M n (MeV) 




CERN-TOV[19] 
ETM [2C, 21] 
JLQCD/TWQCD [||] 


Wilson +0(a) impr. 
Wilson TM 
Neuberger 


0.052-0.072 (M K ) 
0.065-0.1 (F n ) 
0.12 (r ) 


>380 
>265 
>290 


3.2 -3.6 
3.2-3.6 
2.9 



Table 1: Simulation parameters of recent Nf = 2 computations. 



Collaboration Dirac operator 



a (fm) 


(MeV) 




0.06-0.18EF*] 


>240 


4 


0.13 [r ] 


>290 


3.7 


0.11 [M n ] 


>330 


4.6 


0.09 [Md 


>160 


2.3 



MILC Q23, 24J] Staggered 



NPLQCD [ 25 ] Domain Wall + Staggered 



RBC/UKQCD [ 26 ] Domain Wall 



PACS-CS [g7|] Wilson + 0{a) impr. 



Table 2: Simulation parameters of recent Nf = 2+1 computations, ro is the so-called Sommer scale 

adopted, the lattice spacing (in parenthesis the quantity that has been used to fix the scale), the 
pseudoscalar meson masses and the volume in terms of the quantity M n L. 

One of the main outcomes of these studies is that M\ is nearly a linear function of m, as predicted 
by the GMOR relation, up to relatively large masses (of the order of m s /2). A further common 
conclusion is that for M % < 450 MeV lattice results seem compatible with prediction of the chiral 
effective theory at NLO. CERN-TOV and JLQCD pointed out however that for F n NNLO effects 
may be significant and this could affect the systematic uncertainty on the determination of I4. 

Several collaborations are performing lattice calculations with 2+1 dynamical flavours, i.e. 
with two light (degenerate) quarks and one heavier quark (which is usually fixed at the physical 
strange mass m s ). In analogy to the Nf = 2 case, we report the relevant simulation parameters in 
Table || 2 . The MILC data have been analysed using SU(3) x SU(3) (partially quenched) rooted 
staggered chiral perturbation theory including analytic NNLO and NNNLO terms. The NPLQCD 
collaboration adopted Domain Wall valence quarks and staggered sea quarks, and fitted the data 
using continuum chiral perturbation theory at NLO and (partial) NNLO. RBC/UKQCD adopted 
a partially quenched setup, while PACS-CS used only unitary points. These two collaborations 
matched their results for pseudoscalar masses and decay constants with SU(3) x SU(3) effective 
theory at NLO and observed poor convergence of the perturbative series around the physical strange 
quark mass. On the other hand, the NLO SU(2) x SU(2) analysis seems to yield more reliable 
results. 

The results obtained for the LECs are all summarized in the Tables [|, |], |5|, ]6| 

3. Other approaches to determine LECs on the lattice 

3.1 LECs from £-regime simulations 

On a finite volume V = L 3 T with L S> 1 /Aqcd> different chiral regimes can be distinguished. 



2 Another collaboration has recently published results for hadron masses obtained with 2+1 lattice simulations 
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Approaching the chiral limit by keeping M K L ^> 1 (like already mentioned in Section defines the 
so-called /^-regime, where finite-volume effects are exponentially suppressed, while mass-effects 
are dominant. Alternatively, one can approach the chiral limit by keeping pt = mLV < 0(1); this 
defines the so-called £ -regime [|8[ |9|], where the Compton wavelength associated to the Nambu- 
Goldstone bosons is larger than the linear extent of the box, M n L < 1. In this case the power 
counting is reorganised such that mass effects are suppressed, while finite-volume effects are en- 
hanced and become polynomial in L , . One of the consequences of the rearrangements is that, 
at a given order in the effective theory, less LECs appear with respect to the /^-regime, and the 
predictions are less "contaminated" by higher order unknown couplings. For instance, the NLO 
predictions of given correlation functions contain only the LO couplings: F, L (FqXo) can then be 
extracted by matching the finite-size scaling of correlators computed on the lattice with the predic- 
tions of the effective theory. 

As an example, the pseudoscalar and axial correlators at NLO can be written as (for Nf = 2) [29] 



C,(f) = ^ f d*x(P{x)P{0)) =L 7 
C A (t) = f d 3 x{A Q (x)A (0)) = 



ap + j2L? bph \T 



V 



aA + FH? A 1 Vr 



(3.1) 
(3-2) 



where h\ (t) = 2 



i) 2 



, and cip,bp,aA,bA are dimensionless functions of }X, L, T . 



Furthermore, in the £ -regime topology plays a relevant role []30p : observables may be defined at 
fixed value of the topological charge, and the dependence on this charge should be also reproduced 
by the chiral effective theory. 



This strategy has been applied recently in two studies. JLQCD [ ]31[ ] computed mesonic correla- 
tion functions at fixed topology in the £ -regime using the Neuberger Dirac operator with Nf = 2 



dynamical quarks. A. Hasenfratz and collaborators [32] adopted improved Wilson fermions com- 
bined with the reweighting technique. The results obtained for F and L are reported in the summary 
Table |3] (8,9). The two-point functions computed in [ [32| ] have been recently reanalysed in [33] in- 
cluding the leading 0(a 2 ) correction inferred via the so-called Wilson chiral perturbation theory 
applied to the £-regime [ |33| , |34|], finding small corrections with respect to the continuum results. 

Two-point mesonic functions have been computed in the effective theory also in the case of 
non-degenerate quark masses, and in particular in the case where some quarks are in the ^-regime 



and others in the £-regime [J35]-|37|]. This may offer new possibilities to extract the LECs from 
lattice data. 



3.2 LECs from eigenvalues distribution 

At LO in the £-expansion, the partition function is equivalent to the one of a chiral Random 
Matrix Theory (RMT) [|3^-41]; it follows that RMT reproduces the same microscopic spectral 
density ps{^,pi) of the chiral effective theory, in terms of two dimensionless variables £ = XLV 
and jit = mLV, where X represents the eigenvalues of the Dirac operator. Moreover, one can extract 
the probability distributions of single eigenvalues [|2|-[n]]; hence it is possible to match the QCD 
low-lying spectrum of the Dirac operator (A / t) ( 2 CD with the expectation values (^k) RMT (where 
k = 1,2.. labels the eigenvalues) in order to extract the chiral condensate E. This method has been 
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applied by many authors; the results obtained by the JLQCD/TWQCD collaboration [15] an by 
DeGrand and collaborators [46] are reported in Table ||] (10, 11). 

This framework can be extended such that the spectrum of the Dirac operator is sensitive also 
to the pseudoscalar decay constant F at LO in the effective theory [07|]. This strategy has been 
adopted by DeGrand and collaborators [48], yielding the results reported in Table |3| (12). 

In another recent work [^] it is shown that the quark condensate can be extracted from suitable 
(renormalizable) spectral observables defined in the /^-regime, for instance the number of Dirac 
operator modes contained in a given interval. The method has been used on Wilson lattice QCD 
and the result for the condensate is also given in Table |3| (6). In correlation to this, an interesting 
development from the side of the effective theory has been investigated in [|50|], where the quark 
condensate and the spectral density of the Dirac operator have been computed by using a technique 
which is able to smoothly connect p- and £ -regimes. 



3.3 LECs from form factors 

The pion electromagnetic form factor, defined as 

(n+ip'W^n+ip)) = (p + p'),Fr(q 2 ) (3.3) 

with q 2 = —Q 2 = (p — p') 2 , can be matched with the predictions of the chiral effective theory. In 
particular, at NLO it turns out to be related to the LEC l 6 for the SU(2) x SU(2) theory [0] and 



to Lg in the SU(3) x SU(3) case [51]. The pion electromagnetic form factor has been computed 



on the lattice and compared with chiral effective theory by the RBC/UKQCD collaboration [52] 
and by the ETM collaboration [|53|], using the lattice parameters already given in Tables [jj |[ The 
latter performed a NNLO chiral fit, which -by using as input the experimental value of the scalar 
pion radius- allowed the extraction of F, E, l\ , h , h , U , h ■ We report some of the results in Tables [|||. 
Moreover, an exploratory study of the scalar form factor has been presented by the JLQCD/TWQCD 
collaboration at the last lattice conference []54]]. 



4. Summary and Conclusions 

In Tables H||, | we summarize the results for the low-energy couplings obtained through 
lattice simulations. 

In particular, Table ^collects the results for the LO LECs, both from /^-regime and £ -regime 
studies. The first error is statistical, while the following uncertainties (if present) are systematic. 
Although errors are still quite sizeable, the general agreement among the different determinations 
is reassuring. PACS-CS results for the condensate £ point towards larger values with respect to 
other collaboration: notice however that in that case perturbative renormalization has been used 
and systematic errors may be important. Table | shows the ratios F /F$, L/Lo; large deviations 
from 1 would indicate strong violations of the OZI rule. More precise determinations are needed 
for these observables in order to draw a definitive conclusion. 

Table || summarizes the results for the SU(2) x SU(2) NLO constants 1$, I4, For some of the 
lattice simulations with Nf = 2 + 1, we report both values obtained through a direct SU(2) x SU(2) 
fit and by converting SU(3) x SU(3) into SU(2) x SU(2) LECs by means of chiral perturbation 
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Collaboration 



N f 



F (MeV) 



Fq (MeV) Z 1 / 3 (MeV) 



,1/3 



(MeV) 





(1)ETM@|J 


2 


86.03(5) 


267(2X9X4) 




(2) JLQCD/TWQCD [||] 


2 


79.0(2.5)(0.7) (+^) 


235.7(5.0X2.0) (+oo 7 ) 




(3) MILC @ 


2+1 




278(1)(±|)(5) 


242(9) (+ 5 7 ) (4) 


(4) RBC/UKQCD [0] 


2+1 


81.2(2.9X5.7) 


255(8)(8)(13) 




(5) PACS-CS 10] 


2+1 


90.3(3.6) 


83.8(6.4) 309(7) 


290(15) 


(6) Giusti, Liischer[fJ§] 


2 




276(3)(4)(5) 




(7) ETM Q 


2 


86.6(4)(7) 


264(2X5) 




(8) JLQCD||l]| 


2 


87.3(5.6) 


239.8(4.0) 




(8) A. Hasenfratz ef a/[[32|] 


2 


90(4) 


248(6) 




(10) JLQCDI0] 


2 




251(7)(11) 




(ll)DeGrandef al^ 


2 




282(10) 




(12) DeGrand, SchaeferQ 


2 


84(5) 


234(4) 





Colangelo, Dtirr [R8J 



phen. 86.2(5) 



Table 3: Summary of lattice results for LO couplings. The condensate E is evaluated in the MS scheme at 
the scale ji —2 GeV. The upper part (1-7) collects results obtained by matching lattice QCD with the chiral 
effective theory in the ^-regime, while the bottom part (8-12) summarizes results obtained in the e-regime. 
The last row is the estimation from a phenomenological analysis. 



theory. Also in this case the good agreement is encouraging. Lattice data for /3, 1 4 agree also 
with phenomenological estimates ^ |5|]. The phenomenological impact of the LECs I3, 14 on the 
s-wave pion scattering length a^{I = 0,2) has been recently discussed by H. Leutwyler [|56|]. 

Finally, Table |6] shows the lattice results for SU(3) x SU(3) NLO constants L4, L5 and the 
combinations (2L^ — L4), (2L§ — L5), at the scale M p = 770 MeV. These values can be compared 



with phenomenological estimates reported in [57]. In this case the situation is less clear, since 
it has been pointed out that SU(3) x SU(3) chiral fits may not be appropriate to describe lattice 
results obtained with Nf = 2 + 1. Notice that in order to extract SU(3) x SU(3) couplings, lattice 
simulations with Nf = 3 light flavours would be more appropriate. 

To conclude, important progress occurred in lattice QCD simulation in the past years. Un- 
quenched computations are reaching pion masses as light as 200 MeV, lattice spacings around 
a ~ 0.07 fm and spatial extends larger than 4 fm, moving towards the "physical point". This will 
allow a precise matching with the chiral effective theory and a deep understanding of low-energy 
properties of strong interactions. By decreasing systematic uncertainties, systematic errors become 
a very important issue and have to be carefully studied, namely by controlling the continuum ex- 
trapolation, the finite-volume effects, the renormalisation and the uncertainty coming from higher 
orders in the chiral effective theory. LECs obtained with independent methods and different dis- 
cretizations tend to point at uniform results within the errors, giving encouraging perspectives for 
the future. 
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Collaboration 



F/F 



L/E 



MILC [24] 



1.15(5) (t?) 1.52(17) (Tj^) 
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u ■ 10 3 



L 5 ■ 10 3 
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1/ 
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Table 6: NLO SU(3) couplings at the scale M p = 770 MeV extracted from Nf = 2+1 lattice simulations. 
In the last line the phenomenological estimate is reported. 
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